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We present an analytic method based on the Hadamard-WKB expansion to cal- 
culate the self-force for a particle with scalar charge that undergoes radial infall in 
a Schwarzschild spacetime after being held at rest until a time t = 0. Our result is 
valid in the case of short duration from the start. It is possible to use the Hadamard- 
WKB expansion in this case because the value of the integral of the retarded Green's 
function over the particle's entire past trajectory can be expressed in terms of two 
integrals over the time period that the particle has been falling. This analytic re- 
sult is expected to be useful as a check for numerical prescriptions including those 
involving mode sum regularization and for any other analytical approximations to 
self-force calculations. 

PACS numbers: 04.25.Nx, 04.30.Db, 04.70.Bw 



I. INTRODUCTION 

This latest version has a new appendix C containing the Erratum in Ref. [41] to the 
version published in Physical Review D [40]. In this version we have made the corrections 
indicated therein. 

Computation of the self-force for a point particle in orbit around a black hole is a topic of 
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active research today 

HAS 

prompted by the preparation of gravitational wave detectors 
such as LISA which are capable of detecting gravitational waves emitted when a compact 
object falls into a supermassive black hole [J, ]5[. 

An exact expression for the self-force in a black hole spacetime has been obtained only for 
the cases of a scalar or electric charge held at rest in a Schwarzschild spacetime 
and an electric charge held at rest on the symmetry axis of a Kerr spacetime [10l . 111 ]. 
Approximate analytical expressions have been obtained using Green's functions for scalar, 
electric, and gravitational charges in various weak field limits 
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calculations have involved 



such as a static charge [17 



;he use of mode sum techniques |15|, ll6] in cases of hi 



.3, El" 1 



'Q. 



lid 



Most other 
i symmetry 



Il9( , radial infall 20|, |2l| , a circular orbit [22], [23 
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25|, or 



a slightly elliptical orbit 241 ]. 

Since the mode-sum regularization procedure has been developed extensively, aiming at 
practical calculations, it is important to find independent ways to check its accuracy and 
reliability. This of course also applies to any other method, numerical or analytical that 
might be developed in the future. To date the checks on various mode-sum regularization 
procedures which we are aware of fall into three categories: 1) comparisons with exact ana- 
lytical results, 2) comparisons with analytical approximations, 3) comparisons of the results 
of one mode-sum regularization technique with those of another. Most of the comparisons 
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271 ] . Comparisons with 



that have been made fall into the third category [21, 
exact analytical results are of course the most reliable but also the most limited. They have 
only been done in the case of static scalar and electric charges held at rest in a Schwarzschild 

n 

spacetime [17j . Comparison with an analytic approximation in the weak field limit has been 
made for a static scalar charge at rest in an axisymmetric spacetime In this paper we 
present results which can be used as a check in the second category. However, unlike most 
previous analytical approximations to the self-force, ours is valid in the strong field as well as 
the weak field limit. Specifically, we consider the case of a particle with scalar charge which 
is held at rest until a time t — and subsequently falls radially towards a Schwarzschild 
black hole. 

In a previous paper 28] (which we will refer to as Paper I) we introduced an approximation 



for the computation of the self-force of a scalar char ge in a Schwarzschild spacetime using the 



Hadamard-WKB expansion [2i 
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31 



32. 



33, 



34 



351 ] which can account for the radiation 



reaction effects from the charge's recent past history. This is an expansion for the tail 
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term part of the retarded Green's function and is valid in a spacetime region close to the 
charge. A similar approach has been taken to compute part of the gravitational self-force 
in an arbitrary spacetime by Anderson, Flanagan, and Ottewill [36]. In Paper I an eighth 
order WKB expansion was used to compute the nonvanishing part of the tail term to order 
(x — x') 6 . Since then we have obtained results to order (x — x') u using a 16th order WKB 
expansion. 



Recently Anderson and Wiseman [37| investigated the convergence of the Hadamard- 
WKB expansion for the self-force. They posited that convergence is obtained as long as the 
point separation in proper time is small enough so that the null geodesies emanating from 
the point at the earlier time have not had time to re-intersect the particle's path. We know 
the absence of caustics is perhaps the least stringent criterion for the validity of a quasilocal 
expansion (e.g., the domain of validity of Riemann normal coordinates). After computing 
this proper time separation for the case of a particle orbiting a black hole in a circular path, 
they use the result of Paper I [Hadamard-WKB expansion to order (x — x') 6 } to evaluate 
the convergence of the series expansion for the self-force. 

Anderson and Wiseman go on to show that the primary constituent of the self-force 
comes from the tail term of the Green's function when the points are too widely separated 
for this expansion to converge. This would seem to rule out the possibility of using the 
Hadamard-WKB expansion by itself to compute the self-force, which is not surprising, as 
the quasilocal expansion has this intrinsic limitation to begin with. 

However, we want to point out that there is at least one case in which it is possible to 
determine the entire self-force using the Hadamard-WKB expansion. If the particle is held 
at rest until a time t = after which it falls radially towards the black hole, then the entire 
self-force can be determined using the Hadamard-WKB expansion so long as the particle 
has not fallen too far from its starting point. While the range over which the particle falls 
must be small, the results can serve as a useful check on other methods of computing the 
self-force. 

In Section[TT]a brief review is given of the general method we use to compute the self-force. 
In Section IIHI the self- force is computed for the case in which a particle is held at rest until 
time t = and subsequently falls towards the event horizon of a Schwarzschild black hole. 
We derive a closed expression for the self-force in terms of integrals over the tail part of the 
retarded Green's function when the points are split by only a small amount. The splitting 
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is between the point where the particle is at a time t > and various other points inside 
the past lightcone of this point. We then use the Hadamard-WKB expansion described in 
Paper I to derive power series expressions for the nonzero components of the self-force. In 
Section [IV] the series are evaluated for a specific example. The results are given in Table 
I and Figure I. Some details of the calculations are given in Appendices A and B and the 
coefficients derived for the series are also given in Appendix B. 



II. REVIEW OF METHOD 



Consider a massless scalar field with a source term consisting of a point particle of scalar 

tion is 3, i 1 



charge q. The wave equati 



□$(x) 
p(x) 



-p(x) 



(1) 



dr 



where r is the proper time and z(t) is the trajectory of the particle. A formal solution to 
this equation can be obtained using the retarded Green's function: 

/oo 
G R (x,z(T'))dr' . (2) 
-oo 

The self-force is given formally by 

U(r) = q[V^(x)l =z(T) . (3) 

This expression is divergent and must be regularized. Quinn ^ has shown that the regu- 
larized expression for the self-force can be split into a local term plus a finite integral over 



the gradient of the retarded Green's function. The latter is often cal 
The Hadamard expansion for the retarded Green's function is 



ed the "tail term r 



32 



33 



G R (x,x') = 6(x,x') { U ( X,X 'h [(r(x,x')] - V ^ X '^ 9[-a(x,x')] 



Art 



57T 



(4) 



1 Note that our conventions differ from those of 0, S] • One can obtain those of Ref . Q by letting $ — * 
$/\/47r, q — ► V47T9, and Gr — ► GrKAk). Those of Ref. 0] can be obtained by letting $ — * — $/\/47r and 



q — > V Airq. 

2 In Eq. (1.1) of Paper I this equation was written with an incorrect factor of o(x,x') in the denominator 
of the first term. The relationships between our definitions of u and v and those of Ref. [2] were also 
given incorrectly. Taking into account the previously mentioned difference in the definition of the retarded 
Green's function, the correct relationship is that the definition of u in Eq. (4) is equivalent to that in Ref. 
[2] while the definition of v differs from that of Ref. [2] by a factor of 1/2. No other equations in Paper I 
are affected and the results remain unchanged. 
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Here 0[— a(x, x')] is defined to be one if the point x' is inside the light cone of the point x and 
zero otherwise, while 6(x, x') is defined to be one if the point x resides in the future of a space- 
like hypersurface involving the point x' and zero otherwise. The quantity a(x, x') is equal 
to one-half the square of the proper distance between x and x' along the shortest geodesic 
connecting them. The function v(x,x') contributes to the tail part of the self-force [2|. It 
obeys the equation 

O x v(x,a/) = (5) 
and is symmetric under the exchange of the two points, i.e. v(x, x') = v(x', x). In a general 



spacetime 



30 



31| 

v(x,x) = —-R(x) . (6) 

The tail term for the self-force obtained from Eqs. (CD) and ([2]) can be written in the form 
(see, e.g., [22]) 

+52 /I(^ zMi L (T / T '' (7) 

It is necessary that tq be chosen so that the Hadamard-WKB expansion for v is valid 
throughout the region of integration of the first integral in Eq. ([7]). 

An expansion for v (x, x') in Schwarzschild spacetime was found in Paper I using a WKB 
expansion for the Euclidean Green's function. The WKB expansion is obtained via an 
iteration procedure and is increased by two orders upon each iteration. From a WKB 
expansion of order (2N) one can obtain an expansion for v(x,x') that includes terms up to 
order (x — x') 2N ~ 2 . For the metric 

ds 2 = -(1 - 2M/r)dt 2 + dr 2 /(l - 2M/r) + rW , (8) 

the expansion for v is of the form 

oo 

v(x,x') = v ijk (r)(t-t') 2i (cosi-iy(r-r') k (9) 

i,j,k=0 

with 

cos 7 = cos 9 cos & + sin 9 sin 9' cos(0 — 4>') ■ (10) 

We have recently found that it is possible to extend this expansion to order {x— x') 2Ar_1 by 
using the fact that v(x, x') is a symmetric function. If one takes {2N — 1) partial derivatives 
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of v(x,x r ) with respect to some combination of the coordinates x a and sets x' = x then by 
symmetry it must be true that taking the same combination of derivatives with respect to 
x' a and setting x' = x results in an equivalent expression. Using the expansion ([9j) results in 
the coefficients of terms of odd powers of (x — x') being expressed in terms of the coefficients 
of terms of smaller even powers. The coefficients of terms of even powers of (x — x') cannot 
be obtained in this way and must be obtained from the WKB expansion. 



III. COMPUTATION OF THE SELF-FORCE 



We now proceed to compute the self-force for a particle of mass m and scalar charge q 
falling radially from rest. We assume that the particle has been held fixed at r = r , 9 = 6 , 
and cf) = <po from t = — oo to t — which corresponds to proper time r = 0. At time t = 
the particle is released, with no initial velocity, and subsequently falls radially towards a 
Schwarzschild black hole of mass M. We shall denote the spacetime point at which it is 
released by yo- Our objective is to calculate the self-force exerted on the particle at a proper 
time r > when the particle is at the point y = (t, r, 6q, 4>q). 

The first term in Eq. ([7]) can be computed using the Hadamard-WKB expansion. The 
results of this computation are discussed below. The second term is more difficult. We next 
show that it too can be computed using the Hadamard-WKB expansion due to the fact that 
the charge is stationary until it begins falling at time t = 0. 

First consider the problem of computing the field & s t a tic at the point y due to a completely 
static charge q at the position r = r ,6 = 6 ,<f) = <ft . The solution to this problem in a 
Schwarzschild spacetime has been given by Wiseman 8|. For the above locations of the 
charge and field point, with a metric of the form ([8]), and taking into account a difference 
in conventions, it can be written as 



1 / 2M 1 

®stauc(y) = r ?Wl • ll 

On the other hand one can use Eq. ([2]), to write the field due to the static charge as 

/oo 
G R [t,r;t'(r'),ro}dr' . (12) 
-oo 

Here and for the rest of this section we suppress the dependence of various quantities on the 
angles 8q and 0o- Dividing the integration region in the same way as was done in Eq. (J7j) 
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and using the Hadamard expansion (j3j) one can write this latter equation as 



®static{y) = T~\/l- — / r ; t', r o] r; t', r )] (ft' 
4tt V r Jo 



- iS/ 1 / v[*,r;f,r D ]d* / +g / G R [t,r;t'(T'),r }dr' (13) 

ovr V ?"o Jo J-oo 

In the first two terms the integration variable has been changed from r' to t'. The two theta 
functions in Eq. (jl]) result in an upper limit for the second integral which is equal to the 
retarded time t R which is given in Eq. (1A12j) . The time t is taken to be the time that it 
would take a particle to fall from ro to r assuming that it starts at rest. Then the third 
term on the right in Eq. (fTBl) is the same, except for the gradient and a factor of q, as the 
second term on the right in Eq. <^§. The value of this term can be obtained by computing 
the other three terms in the equation. The term on the left is given in Eq. fill I) and the 
second term on the right can be computed using the Hadamard- WKB expansion. 

To calculate the first term on the right in Eq. ([TBI , we note that the argument of the 
delta function vanishes on the light cone of the point y. Since the charge is static 

a(t,r;t',r ) = a(t,r;t R ,r ) + (^pcr(t,r;t' ,r )j \ v=tR {t'-t R ) + ... 

= a tR (t?-t R ) + ... , (14) 

where the shorthand notation = cr ;At has been used. Then 

5(a(t, r] t R ,r )) = 6[a tR (t'-t R )] = ^ - tR) . (15) 

Wt R \ 

Next one must calculate u(y,y R ) with y and y R = (t R , r , 9 , O ) connected by a null 
radial geodesic. By substituting the Hadamard expansion into the equation satisfied by the 



Green's function it is possible to show that in general 



29. 



30 



3l| 



u(x,x') = A 1/2 (x,x') (16) 
A(x,x') = 



det(-o- ;Atl/ ) 



y/-9(x)y/-g{x') 

Thus what remains is to calculate a tR and a ;Atl / for the two points x = y and x' = y R . 
Although there may be some simple way to reason out the answer, as shown in Appendix 
A, it can be obtained by solving the geodesic equations and integrating the result to obtain 
the proper distance along the geodesic. The result is 

a tR {t,r;t R ,r ) = r - r 
u(t,r;t R ,r ) = 1 . (17) 
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Substituting Eq. (fTTI) into Eq. (fT3|) and computing the integral one finds that 

/ G R [t,r;t'(r'),r }dT' = —Jl / v[t, r; t', r }dt' (18) 



With the definitions 



®s(y) = ^-Jl-— / v[t,r;t',r ] dt' (19a) 



vr V r J 



®f(v) = ^ / v[t,r;t'(r'),r'(T')}dr' (19b) 



Eq. ((Tj) becomes 



*-(*.(*)-*/(*)) 



(20) 



Here the facts that in a Schwarzschild spacetime v(x, x) — and (as shown by our expansion) 
= 0; have been used to interchange the order of integration and differentiation. Note 
that this derivation only works for the time and radial components of the self-force. Because 
of spherical symmetry, the angular components of the self-force for a radial trajectory are 
zero. Finally the subscript "tail" has been dropped because for a geodesic trajectory the 
local part of the self- force is zero in a Schwarzschild spacetime 

As a result of Eq. (120]) . the problem of calculating the self- force reduces to calculating $ s 
and $/. This is an exact result. We now calculate the right hand side of Eq. (J20l using the 
Hadamard-WKB expansion for v(x, x') whose form is given in Eq. ([9]). For radial geodesies, 
cos 7 = 1, so only the coefficients v i0 k(r) contribute. The result for <3> s is 



= i E (aTl) Wr) - (t - t R f^] (r - r f (21) 

i,k=0 



To calculate $/ one can use the geodesic equations (1A7I) to convert the integral (119bj) to 



an integral over the radial coordinate r. One can further solve the geodesic equations to 
obtain the trajectory t(r). After substituting the Hadamard-WKB expansion for v, the 
integral ( 119bl) can be computed numerically. 

An alternative is to expand all relevant quantities in both <3> s and $j in Taylor series about 
r . This allows one to compute the integrals analytically order by order. The derivation is 
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given in more detail in Appendix B. We find 
ft{r) = ^.t ., Jl - — ( 10- 



fr*(r) 



22400 7T r\ 

3q 2 
11200^ 



r 



106M\ fr 
r J 



2M 
r 



2M 
ro 



r — r 
r 



ro 

5/2 



+ 



r — r 



r 



O 



r — r 
r 



7/2 



with r* the Regge- Wheeler coordinate defined by 

V - 2M 
2M 



(22a) 
(22b) 

(23) 



r* = r + 2M log 

It turns out that each subsequent order of the WKB expansion adds another term to the 
series. Using a 16th order WKB expansion we have results for a total of six terms in the 
expansions for both f t and / r *. The coefficients of these terms are displayed in Appendix B. 



IV. SPECIFIC EXAMPLE 



As a specific example we have chosen the case in which a particle begins falling from 
rest at r* = 40M at time t = r — 0, with r* defined in Eq. ( T23l) . Using the series in 
Appendix B, the ft and f r * components of the self-force have been computed. The results 
are displayed in Table I and Figure I. The relative errors in the table are estimated by taking 
the absolute value of the ratio between the last term used in the series expansion fICll) for 
a given component of the self-force and the entire series for that component. 

The expected convergence distance, obtained by measuring how far the particle will fall 
in the time for a null geodesic leaving the location of the particle at time t = to circle 
the black hole and intersect the particle in a new location, is r « 30. 7M. On the other 
hand, an examination of the error shows that more terms than the six we have computed for 
each component of the self-force are necessary for an accurate determination of the self-force 
when r < 33M. 

While the range over which we can compute the self-force is relatively small and the calcu- 
lation is for a particle held at rest that subsequently undergoes radial infall in a Schwarzschild 
spacetime, the merit of this method is it produces reliable results and thus can provide an 
independent check on other (existing and future) prescriptions for the calculation of the 
self-force. It is worth noting that our calculation started with a finite result and therefore 
no regularization was needed. The generalization to the cases where the self-force is due to 
electromagnetism and gravity should be straight-forward. 
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r/M 


W l4 (M 2 /q 2 )f t 


Error 


10 14 (M 2 /g 2 )/ r , 


Error 


34.25 


6.93066 


6.9 x 10~ 6 


-5.7384 


3.4 x 10~ 5 


34.00 


120.72 


4.4 x 10~ 4 


-96.78 


1.4 x 10~ 3 


33.75 


616 


3.4 x 10~ 3 


-520 


8.4 x 10~ 3 


33.50 


2080 


0.013 


-1840 


0.25 


33.25 


5680 


0.031 


-5200 


0.052 


33.00 


13800 


0.059 


-13200 


0.089 


32.75 


30000 


0.10 


-30000 


0.13 


32.50 


66000 


0.14 


-62000 


0.18 



TABLE I: The dependence of the temporal (ft) and radial (f r *) components of the self-force on 
the radial distance r/M are given for a particle undergoing radial infall after being held at rest 
until a time t = 0. The particle's initial location is at r* = 40M which corresponds to r ~ 34.43M. 
The error shown for each case is an estimate of the relative error and is obtained by taking the 
absolute value of the ratio of the last term used in the series for the self-force to the entire series. 
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APPENDIX A: CALCULATION OF u AND a tR 



In this appendix we calculate the quantities u(y, i/r) and (Jt R (y, Vr) with y = (t, r, #0, 4>o) 
and dr = (tn, r , #0, 0o); tn is the retarded time defined in Eq. (1A12I) . The points y and yR 
are separated by a radial null geodesic. From Eqs. ffT4l) and f|T6|) it is clear that 

1/2 

(Ala) 



u(y,y R ) 



-det (—a 



; iiu' 



(x,x')) 



9 i >\ 



(Alb) 



Although the points are separated by a radial null geodesic, it is necessary to assume a 
more general separation before computing derivatives. Once this is done then the specific 
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33.0 



33.5 



34.0 



54.5 



r/M 



FIG. 1: In this figure the two components of the self- force are plotted for the case of a particle 
undergoing radial infall after being held at rest until a time t = 0. The particle's initial location is 
at r* = 40M which corresponds to r « 34.43M. The dashed curve corresponds to ft and the solid 
curve to f r * . 

separation can be taken. Because of spherical symmetry the angular dependence of the 
quantity a(x,x') must be such that it is a function of cos 7 which is defined in Eq. (flOl) . To 
see this assume that the two points have an angular separation of 7. Since any arbitrary 
rotation of the coordinate system does not change this angle, the angular dependence of 
a can only be through 7. Because any 2tt rotation around the origin of the coordinate 
system should render a unchanged, a must be a periodic function of 7 with period 2tt. This 
dependence can be written as a Fourier series in sin 117 and cos 717. The coefficients of the 
sine terms must identically vanish since a(x, x') is a symmetric function in x and x' and 
switching x and x 1 changes the sign of 7. Since cos 717 can be written in terms of powers of 



cos 7 



it is then clear that a = a(t, r; t', r'; cos 7). 
We next compute the angular derivatives of a beginning with 

da 



a.a[x,x 



[— sin 9 cos 9' + cos 9 sin 9' cos(</> — <$')) 



(A2) 



s <9cos7 / 

Note that in general cr ;At (x, x') is a vector at x and a scalar at x' . Next we can proceed to 
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take a second derivative. 

a-00i(x, x') = | ttt \~o ) ( — sin 9 cos 9' -\- cos sin 9' cos(0 — 0')) 2 

\a(cos 7 )^/ 

+ ( — — — J (sin sin 0' + cos 5 cos 0' cos(0 - <p')) (A3) 
\ac0s7 / 

For the special case of a radial geodesic, 9' = 9, <p' = <p, and thus cos 7 = 1. As a result 

faw) BM ~ =1 =(J^—) ■ (A4a) 

7 V^cos7; cos7=1 

Similarly it is easy to show that 

(a. ul ) , = sin 2 9 ( — da - 1 (A4b) 

and 

( cr ;^') C o S7 =i = ( (T ;er')cos 7 =i = (<W) C o S7 =i = ( A4c ) 

(°";^')co S7 =l = ( a ;^') C os 7 =i = (°";0e')cos7=i = • ( A4d ) 

Along with the other combinations of radial and time derivatives it is still necessary to 
calculate ( tJ^ 3 - > ) . If we choose coordinates such that 9 = 9' = tt/2 then from Eq. 

V ''J cos 7 =1 

(EibD 

( 0/ - V J = ( <T \<w)<i/=& ■ (A5) 

To calculate the rest of the derivatives that comprise a. ^ it is useful to write a in terms 
of the proper time r and to assume that the points both lie in the equatorial plane so that 
9 — 9' — tt/2. Since a is one-half the square of the proper distance between the two points 
along the shortest geodesic connecting them we have 

a = -V 2 . (A6) 



The proper time r can be computed by solving the geodesic equations 



r = -y/E / jTjT- (A7a) 

r dr 

t-t' = - «77^ (A7b) 

Jr> f(f)(l-Ef(f)-^y/l 



J I f 2 (1 _ E/(f) _ W))l/2 (A?C) 
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where f(r) = 1 — 2M/r. Note that J = yields a radial geodesic and E = a null one. 

Expanding the integrand of Eq. (IA7al) in powers of E and J, computing the integral, and 
then substituting into Eq. (1A6I) gives 

a = -|V - r) 2 - y F ( t> ~ r ) + ^ 3 ) + °(^ j2 ) ( A8 ) 

with 

F = - dff(f) = r' - r - 2Mln . (A9) 

The values of the constants E and J can be found by solving Eqs. (1A7bl) and 0A7c0 . The 
integrals can be computed by first expanding the integrands in powers of E and J. From 
Eq. (1A7cl) one finds 

J = vzr r U> - + °W ~ ^ + °U ~ ^') 3 ] • ( A1 °) 

Keeping terms in the expansion of the integrand in Eq. (lA7bj) to 0(J 2 ) and 0(E 2 ) and using 
( lATOj) yields 

E = 2 Jf-lM rV U $f Ht'-tn) 2 F 

- ^) 3 ] + 0[(t' - t R )(</> - <p'f] + OU - <P'f] (AH) 

with 

[ r df , , x r'-2M /A N 

t R = t+ = t -(r'-r)-2M\n — . A12 

7r' f{ r ) r-2M 



Substituting into Eq. (IA8I) gives 

a = (f - t fl )(r' - r) - (f - t R f + ± (0 - 0') 2 

+0[(f - t*) 3 ] + 0[(t' - ^)(0 - 0') 2 ] + - 0') 4 ] • (A13) 
Computing the various derivatives of a and then setting t' = t R and </>' = </> gives 

o-* B = (r' - r) (A14) 

and 

0"; it' 



C;t r ' = -1 + 



f(r')(r' — r) 
F 

~~ _1 + Trvn — V 

j(r)(r — rj 
1 1 

(J. rr i = — — — — - — h 



f(r) f(r>) f(r)f(r')(r'-r) 
= -rr' . (A15) 
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Using Eqs. (1A4j) and (1A15j) it is easy to show that 

det[-cr ; ^,(y,y R )) = (rr sin 6 ) 2 . (A16) 
Substituting into Eq. (lAlaj) then gives 

u(y,y R ) = l. (A17) 

This is an exact result. 

APPENDIX B: POWER SERIES EXPANSION FOR THE SELF-FORCE 

In this appendix an expansion for the self-force in powers of r — r is derived. The 
coefficients of the expansion depend on the mass M of the black hole, the radius r at which 
the particle begins falling, the radius r which is its present location, and the charge q. 

To begin, consider $ s in either Eq. 11 19 aft or ( l2Tj) . Its contribution to the self-force is given 
in Eq. ( 1201) . From Eq. ( 1A12I) it is clear that d t = d tR . Because v(x, X\) satisfies the equation 
n x v = and because of the time translation and time reversal invariance of the metric ([8]), 
v is a function of (t — t{) 2 . Making use of this fact one finds that 



d$ s q 2M 



dt 8n V ?"o 



tR Q 

v{t, r; t R , r ) - / r ! *i> r o) dt i 



= f\l-—v(t,r;0,r ). (Bl) 
8tt V r 

Here as in Section (Till we suppress the dependence of v on 9 and 4> . 

Next consider $/ in Eq. ( 119bl) . The geodesic equations ( 1A7I) can be used to change the 



integration variable from the proper time to the coordinate time. In this case J = as the 
geodesic is radial. Since it starts from rest at r it can be seen from the geodesic equations 
that E = (1 - 2M/r )" 1 . The result is 

^ f = ±—^=[ v (t,r;t 1 ,n)(l-— y ]dt 1 . (B2) 



Noting that in a Schwarzschild spacetime v(x,x) = 0, one finds 

, i [<(dv(t,rM,n ) \( l _2 E dt ^ 



Ot 8lT ft 2M J \ &ti J \ T\ 



'■() 
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Since the particle is freely falling, r\ = and one can write 



dv f 2M 
dti V n 

with the result that 



d 



v 1 



dt 



8n 2M 
ro 

_d_ 



+ 



Thus 



ft 



•57T 



i _ Mi J ro dn 



2M 
n 

v[t,r; 0,r ] ^1 - 
w(t,r;ti,ri) ^1 - 

d 



dr 1 d 
dti dri 



2M 
ro 

2M 

n 



v 1 



2M 

n 



(B4) 



dr\ 



(B5) 



2M 



dri . 



(B6) 



The computation of f r * is straightforward. Taking the derivative of Eqs. ( 12TI) and ( 119bl) 
and then using the geodesic Eqs. (1A7I) and (I2"5j) one finds 

^2 



dv iQk 



2M 



( 1 r ) \ 2i + 1 

(r - r ) fe + fc^ fe(r - r ) k ~ l 



dr 

2 / 2M 



[t 2i+1 ~(t- t R ) 2i+1 ] 



87T 



r 



r r 9f(t,r;ti,ri) / n 



2M 



<9r 



-o?ri . 



(B7) 



fro - V r - r l 

The next step is to expand t in powers of r — r and t\ in powers of r — ri. This is done 
using the geodesic equation (1A7bl) . Then making the change of variables 



ro — r 



r 



si 



x 



W Q 



r - n 



r 



2M 
r 



1 - 



2M 



gives 



ro 

Xq 



2s 
w 



w, 



s 

Wo 



(B8) 



(B9) 
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Substituting this and the corresponding expression for t 1 into Eqs. (1B6I) and (1B7|) . expanding 
in powers of s and s\, and computing the integrals gives 

oo 

ft = q 2 J2 a ^ s2n ( B10a ) 

n=3 
oo 

fr* = <7 2 $>2n + l* 2n+1 • (BlOb) 

n=2 

Using a 16 th order WKB expansion for v (x, x') we find 
w (-13 + 53 w 2 ) 



a 6 



22400 7rr 



2 



fl8 = , 7 n,nrT 2 2 (" 703 + 2031 w o 2 + 72 w 4 ) 
470400 7r ?V £0 

fl io = 1 ^oQon^n 2—4 H 04385 + 1485889 W - 1990813 w 4 + 149 7309 w 6 ) 

au = ~„ 17n ™° h—s (-1380534784 + 6209685607 W ~ 9489117981 w 4 

355170816000 7ir x b 

+2956643877 w 6 + 3255643281 w s ) 

a u = — 5 — a (-690105799163 + 3725808468235 w 2 

129282177024000 7rr 2 x 8 V 

-7969067122446 w 4 + 9871663249718 w 6 - 9633792321223 w 8 

+5317042452879 w 10 ) 

a w = — — - (917523531344444 - 5807610696686037 w 2 

131867820564480000 vrr 2 x w v 

+15120303760167717 w 4 - 19548426682161946 w 6 + 7163353443413246 w 8 
+11432898205974399 io 10 - 9956772991427823 w 12 ) (BlOc) 
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and 



w 4 



11200 7rr^ 

b 7 = — — (-147+ 1031 w 2 - 2004 w 4 ) 

470400 7rx r 2 v ; 

b 9 = ^ (725274 - 4840427 wo 2 + 7544064 w 4 - 1384911 w 6 ) 

558835200 nxfrfi v ' 

b n = -J =-5 (-5192394350 + 35645887586 w 2 - 76724155827' 

1598268672000 irx 5 Q r 2 v 

+74055537336 w 6 - 34933610745 u> 8 ) 

6 i3 = ^oK^ooLnn —( 8865041861 °- 652839443586 W + 1715470441205 w 4 

13851661824000 Tra^r^ 

-2048237904519 u? 6 + 938994638717 w 8 + 27641613573 «; 10 ) 

6i 5 = 5-5 (-1082712168000450 + 8690022729803252 u> 2 

5 98900865423360000 Trxg v 

-27039642594514215 w 4 + 43159798466443548 w 6 - 40035939562445204 w 8 
+24271024582691064 w 10 - 8499901220531595 w 12 ) . (BlOd) 
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APPENDIX C: ERRATUM 



In the Physical Review D (PRD) version of this paper [40j there were errors due to a 



factor of 2 error in Paper I [2 
also published in an erratum 
here as wel 







. These errors have been corrected in this arXiv version and 



411 ] . For completeness we list the corrections in that erratum 



In Ref. 



281 ] there was a factor of two error in the calculation of the coefficients of the 



expansion for v (x, x') that has been corrected in Ref. 42]. This expansion, carried to higher 
order, was used in the calculations of Ref. 40]. We correct the resulting errors in Ref. 
here. 

Eqs. (22. a) and (22. b) should be multiplied by a factor of 2 with the result 



fr*{r) 



22400 vrrg 

3q 2 
112007rrg 



r 



?"0 / 



2M 
r 



1 - 



2M 
r 



ro — r 
r 



r 

5/2 



ro — r 



r 



O 



r — r 
ro 



7/2" 



(Cla) 
(Clb) 



The entries in columns 2 and 4 of Table 1 should be multiplied by a factor of 2. The 
revised results are shown in Table 1. The curves in Figure 1 need to be multiplied by a 
factor of 2. The revised figure is shown in Figure 1. 
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r/M 


W l4 (M 2 /q 2 )f t 


Error 


W l4 (M 2 /q 2 )f r * 


Error 


34.25 


6.93066 


6.9 x 10~ 6 


-5.7384 


3.4 x 10~ 5 


34.00 


120.72 


4.4 x 10~ 4 


-96.78 


1.4 x 10~ 3 


33.75 


616 


3.4 x 10~ 3 


-520 


8.4 x 10-3 


33.50 


2080 


0.013 


-1840 


0.25 


33.25 


5680 


0.031 


-5200 


0.052 


33.00 


13800 


0.059 


-13200 


0.089 


32.75 


30000 


0.10 


-30000 


0.13 


32.50 


66000 


0.14 


-62000 


0.18 



TABLE II: The dependence of the temporal (ft) and radial (f r *) components of the self-force on 
the radial distance r/M are given for a particle undergoing radial infall after being held at rest 
until a time t = 0. The particle's initial location is at r* = 40M which corresponds tor w 34.43M. 
The error shown for each case is an estimate of the relative error and is obtained by taking the 
absolute value of the ratio of the last term used in the series for the self-force to the entire series. 

Equations (BlOc) and (BlOd) need to be multiplied by a factor of 2. The result is 



a 6 
a 8 
aio 



a 16 = 



w (-13 + 53 w 2 ) 
22400 7rr 2 

A7 ™° 2 2 ("703 + 2031 W + 72 w 4 ) 
470400 irr z x 2 ' 

— — j (-404385 + 1485889 W - 1990813 w 4 + 1497309 w 6 ) 

155232000 7rr 2 x 

o^ 17nfi1R nnn — (-1380534784 + 6209685607 w 2 - 94 89117981 w 4 

355170816000 irr 2 x b v 

+2956643877 w 6 + 3255643281 w 8 ) 

— — — - (-690105799163 + 3725808468235 w 2 

129282177024000 irr 2 x 8 v 

-7969067122446 w 4 + 9871663249718 w 6 - 9633792321223 w 8 
+5317042452879 w 10 ) 

5 o — ttt (917523531344444 - 5807610696686037 w 2 

131867820564480000 7rr 2 x 10 v 

+15120303760167717 w 4 - 19548426682161946 w 6 + 7163353443413246 w 8 
+11432898205974399 w w - 9956772991427823 w 12 ) (C2a) 
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^ -5.0x10 - 



-1 .Oxio 4 1 



1 .5x1 



34.5 



FIG. 2: In this figure the two components of the self-force are plotted for the case of a particle 
undergoing radial infall after being held at rest until a time t = 0. The particle's initial location is 
at r* = 40M which corresponds to r « 34.43M. The dashed curve corresponds to f t and the solid 
curve to f r * . 



and 

h 
b 7 

h 

fen 

bi3 



bis = 



3x Q w 2 
11200 vrrg 

5 (-147+ 1031 w 2 - 2004 w 4 ) 

470400 Trxorg v ^ 

ccooo * 3-2 (725274 - 4840427 wo 2 + 7544064 u> 4 - 1384911 w 6 ) 

558835200 irxfro v 7 

=-5 (-5192394350 + 35645887586 w 2 - 76724155827 

1598268672000 7tXq Tq 

+74055537336 w 6 - 34933610745 u> 8 ) 

iQ^i^iJ/mnn ~ (§8650418610 - 652839443586 w 2 + 1715470441205 u> 

13851661824000 7TXgrQ v 

-2048237904519 w 6 + 938994638717 w 8 + 27641613573 w 10 ) 

1 5-5 (-1082712168000450 + 8690022729803252 w 2 

98900865423360000 tt x§ v 

-27039642594514215 w 4 + 43159798466443548 w 6 - 40035939562445204 w 8 
+24271024582691064 w 10 - 8499901220531595 w 12 ) . 
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